Preparing pseudo-pure states with controlled-transfer gates 
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The preparation of pseudo-pure states plays a central role in the implementation of quantum 
information processing in high temperature ensemble systems, such as nuclear magnetic resonance. 
Here we describe a simple approach based on controlled-transfer gates which permits pseudo-pure 
states to be prepared efficiently using spatial averaging techniques. 
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Nuclear Magnetic Resonance (NMR) spin systems pro- 
vide a simple approach for demonstrating techniques for 
quantum information processing [l[ . Because of the diffi- 
culty in preparing pure spin states almost all experiments 
have used pseudo-pure states 0, Q j sometimes called ef- 
fective pure states [4], Q . These are mixed states of the 
form 
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where |0) = 00 . . . 0) is the desired ground state and 
1/2" is the maximally mixed state. The effective pu- 
rity 6 of these states is very low if prepared from high- 
temperature thermal systems @, but the pseudo-pure 
state approach remains useful for simple demonstrations, 
and may also find applications in systems prepared at 
spin temperatures which are low enough to substantially 
enhance the ground state population but not low enough 
to generate a true pure state. 

A range of methods for preparing pseudo-pure states 
have been described which can usually be divided 
into temporal averaging approaches, in which the final 
result is averaged over a number of repetitions of the 
experiment with different starting states, and spatial av- 
eraging approaches, in which magnetic field gradients are 
used to average the spin system over a macroscopic sam- 
ple. Spatial averaging has the advantage that only a sin- 
gle experiment is required, but suffers from a number of 
practical disadvantages: in particular it can be difficult 
to design experiments which prepare pseudo-pure states 
with the highest possible purity, and to avoid the gener- 
ation of zero-quantum coherence terms Q , which are not 
averaged by field gradients. Here we describe a simple 
spatial averaging approach which prepares pseudo-pure 
states with the highest possible purity. 



I. CONTROLLED-TRANSFER GATES 

The thermal state of an NMR spin system is diagonal 
in the computational basis with a pattern of populations 
determined by Boltzmann factors. Although details vary 
greatly, all methods for preparing pseudo pure states fun- 
damentally work by averaging the populations of states 
other than the ground state, ideally leaving the ground 
state untouched @, [9| . 

Averaging of this kind can be achieved using a 
controlled-transfer gate, which comprises a controlled- 
rotation gate with a rotation angle of 9, followed by the 
application of a magnetic field crush gradient which re- 
moves all off-diagonal terms (if the system's density ma- 
trix is diagonal before the gate zero-quantum terms will 
not be generated). The controlled rotation occurs around 
some axis in the xy plane; the choice of axis is unim- 
portant as the off-diagonal terms, which depend on the 
choice of axis, are removed by the field gradient. For 
dcfiniteness the rotation can be thought of as occurring 
around the x-axis, so that a rotation with = 7r corre- 
sponds to a controlled-NOT gate. In the remainder of this 
paper we use the network symbol 



(2) 



to denote a controlled-transfer gate, with the crush gra- 
dient being implicit. Related ideas have been explored 
using line-selective pulses [13, [H| • 

The action of a controlled-transfer gate on a general 
population state of a two-qubit system is 
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c + d + (c — d) cos ( 
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II. HETERONUCLEAR SPIN SYSTEMS 



so that the populations of two states are partially aver- 
aged, with the difference in the two populations being 
scaled down by cos 0. Interchanging the roles of the con- 
trol and target qubits serves to average populations b and 
d, and so the network 
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will mix the populations b, c and d, while leaving a, the 
population of the ground state, untouched. The final 
state will depend on the choice of rotation angles and the 
initial state; to generate a pseudo-pure state the angles 
must be set to 



i\ = arccos 



2b - (c + d) 
3(c - d) 



(7) 



where the gates in Eq.[6]are applied from left to right, and 
the initial density matrix, Eq. [3J is defined with the state 
of the lowest qubit varying most rapidly. If the order of 
the two gates is interchanged to give the network 
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then the values of b and c must be interchanged in Eq. [7] 
For a thermal state of a homonuclear spin system, so 
that the polarizations of the two nuclei are the same, the 
traceless part of the density matrix, sometimes called the 
deviation density matrix [l2j , is proportional to <r\ + a\ . 
In this case d — —a and b = c = so 0i = arccos(l/3) « 
70.5°. As the desired state 1 00) has the largest popula- 
tion of any spin state at thermal equilibrium, and the 
controlled-transfer gates leave this population term un- 
touched, this approach will automatically generate the 
desired pseudo-pure state with the largest effective pu- 
rity which is possible by any averaging process. Meth- 
ods based on exhaustive temporal averaging Q will also 
reach this limit, but the standard methods based on spa- 
tial averaging will not achieve this. 

Note that as controlled-transfer gates only act to trans- 
fer population terms, and these transfers can be calcu- 
lated analytically, there is no need to perform a full sim- 
ulation of the whole density matrix. Although this point 
is not particularly important in two-qubit systems, it be- 
comes important in larger systems as the size of the den- 
sity matrix rises very rapidly with the number of qubits. 
This makes a full simulation impractical, while calculat- 
ing the effects of controlled-transfer gates remains rela- 
tively tractable. 



The method of controlled-transfer gates makes no as- 
sumptions about the initial density matrix beyond the 
fact that it is diagonal, and so can be used with any 
initial population state. In heteronuclear spin systems, 
or spin-systems not starting at thermal equilibrium, it 
is of course necessary to use the more general formula 
in Eq. [7] In some cases only one of the two networks, 
Eqns. |6] and [U can be used, with the other network re- 
quiring an impossible value for 9\ . For heteronuclear spin 
systems starting at thermal equilibrium a solution can al- 
ways be found if the target of the first controlled-transfer 
gate is the spin with the larger polarization; the angle of 
the first gate then rises from 0i = arccos(l/3) « 70.5° to 
0i = 7T — arccos(l/3) « 109.5° as the polarization ratio 
is increased from unity towards infinity. 

In heteronuclear spin systems the controlled-transfer 
approach gate approach can give even larger improve- 
ments in effective purity compared with the conventional 
spatial averaging techniques, firstly because some con- 
ventional spatial averaging methods 13] trade simplicity 
for efficiency, and secondly because such methods usually 
begin with an initial sequence that equalizes the polar- 
ization of the two spins. Although an efficient method 
for doing this in a two spin system is known [l]|, the 
simplest general method in larger spin systems is just 
to reduce the polarization of the more highly polarized 
spins to that of the least polarized spins, inevitably wast- 
ing polarization in the process. 

A special case occurs when the polarization of one 
of the spins is twice that of the other, as in this 
case a pseudo-pure state can be prepared with a single 
controlled-transfer gate with an angle of 7r/2 applied to 
the more polarized spin. This individual result has been 
known for may years, and has applied to homonuclear 
spin systems where the polarization of one spin has been 
reduced [3J and to heteronuclear spin systems where the 
excess polarization of the more polarized spin has been 
reduced to a factor of two [H, [la|; this is not, however, 
a true example of the controlled-transfer gate approach 
as it does not prepare states with the highest possible 
purity. 



III. IMPLEMENTATIONS 

Controlled-rotation gates can in principle be imple- 
mented as a pair of Hadamard gates applied to the 
target spin separated by a controlled-phase-shift gate. 
Hadamard gates can be replaced by pseudo-Hadamard 
gates as usual p}, while controlled-phase-shift gates cor- 
respond to a combination of evolution under the spin- 
spin coupling for some appropriate time, with undesired 
interactions being suppressed by spin echoes (l6j , and lo- 
cal z rotations. If desired, extraneous z rotations and 
cou plin gs can be absorbed into the spin reference frames 
[l7l Il8|| , and any such rotations remaining when the gra- 
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FIG. 1: Experimental demonstration of pseudo-pure states prepared using a conventional approach and our new approach. The 
left hand spectra (a) show thermal states, while the central spectra (b) show pseudo-pure states prepared using a conventional 
approach [la ], and the right hand spectra (c) show pseudo-pure states prepared using our new approach. The small signal 
visible at the centre of the 1 H spectrum arises from an impurity (unlabelled formate). As absolute intensities are not meaningful 
in NMR spectra, the 1 H and 13 C spectra in (c) were normalised to have the same height and all other spectra plotted on 
the same vertical scale. 



dient pulse is applied can simply be dropped, potentially 
significantly simplifying the implementation. Alterna- 
tively controlled-rotation gates can be implemented di- 
rectly using methods such as GRAPE sequences [13, [20j . 

A controlled-transfer gate comprises a controlled- 
rotation gate followed by a crush gradient pulse to remove 
all off-diagonal terms from the density matrix. Imple- 
mentation of a single crush gradient is straightforward, 
but with two or more controlled-transfer gates it is neces- 
sary to consider the possibility of echoes from the interac- 
tion between two successive gradients separated by pulses 
Q • These can be minimised by using gradients with dif- 
ferent integrated strengths or along different axes, and 
will also be suppressed by decoherence and diffusion [2l| . 

Alternatively the gradients can be replaced by tem- 
poral averaging, as applying a crush gradient is equiva- 
lent to averaging spectra acquired with and without the 
application of Z gates (180° rotations); when used with 
controlled-transfer gates it is only necessary to apply Z 
gates to the target qubit. As before the Z gates can be 
implemented as frame rotations, and so this approach 
can be implemented by phase cycling [7] , without chang- 
ing the underlying pulse sequence. This approach allows 
controlled-transfer gates to be used in experimental sys- 
tems where spatial averaging is not available, but as the 
number of experiments required doubles with every gra- 
dient pulse replaced by Z gates the method is only practi- 



cal for networks with small numbers of controlled-transfer 
gates. 

Experiments were implemented using the two-qubit 
heteronuclear spin system ( 1 H and 13 C) provided by dis- 
solving 13 C labelled sodium formate in D2O [2^]. All ex- 
periments were performed at 20° C, on a Varian IN OVA 
spectrometer with a nominal *H frequency of 600 MHz. 
Both nuclear spins were placed on resonance in their re- 
spective rotating frames, so the spin Hamiltonian con- 
tains only a scalar coupling term, with J = 194.7 Hz. 
We chose to implement our controlled gates using se- 
quences of pulses and delays, but all z-rotations (both 
desired rotations and extraneous rotations arising from 
chemical shift evolution) were absorbed into the reference 
frame. The observed relaxation times were Ti = 5.5 s 
and T 2 = 0.61s for the X H qubit and Ti = 14.7 s and 
T2 = 0.35 s for the 13 C qubit. All pulses used were naive 
rotations rather than composite pulses such as BB1 [23| . 

The results, depicted in Fig. [TJ show that this method 
can indeed be used to prepare pseudo-pure states. Both 
1 H and 13 C spectra are shown after applying a 90° exci- 
tation pulse to the spin subsequently observed; no signal 
is seen if no excitation pulse is applied or the excitation 
pulse is applied to the other spin (data not shown). For 
a perfect pseudo-pure state each spectrum should show 
a single absorptive line on the right hand component 
of the doublet, with no signal in the other component, 
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and the height of this line indicates the efficiency with 
which the pseudo-pure state has been prepared. Any off- 
diagonal elements in the density matrix would be visible 
in spectra acquired without excitation pulses, or with 
these pulses applied to the wrong spin. The absence of 
signal in these spectra (data not shown) indicates that 
the final crush gradient has been effective in removing 
off-diagonal terms. Any deviations in the diagonal terms 
from the desired pattern for a pseudo-pure state, that 
is any imbalance in the three population terms suppos- 
edly equalised by the controlled-transfer gates, is visible 
as non-zero intensity on the left hand components of the 
doublets. Pleasingly the error terms are visibly smaller 
with our new method, and the expected increase in signal 
size over the conventionally prepared state is indeed 
seen. 



IV. ASYMPTOTIC APPROACHES 

Although the network described above, Eq. [51 is rela- 
tively simple, it will in general require controlled-transfer 
gates with two different angles, and the initial state must 
be accurately known to determine these angles. It is use- 
ful to consider what can be achieved with simpler net- 
works such that #i = 62 = 0. We will initially confine our 
discussion to homonuclear systems beginning at thermal 
equilibrium, and consider the simplest case, 9 = tt/2. 

With this restriction it is not possible to prepare a 
pseudo-pure state perfectly with a finite number of gates, 
and so it is necessary to have some criterion to measure 
the quality of approximate pseudo-pure states. Conven- 
tional fidelity definitions are not appropriate, as these 
depend only on the population of the ground state [1] 
and are not affected by controlled-transfer gates. Instead 
we use the square root of the x 2 difference between the 
diagonal density matrix p and that of the desired pseudo- 
pure state po , divided by the corresponding value for the 
initial thermal state pi 



e(p) 



tr[(pi - po) 2 



\\P~ Po\\f 
\\pi - Po\\f' 



(9) 



where ||M||f is the Frobenius norm of M, and seek the 
state with the lowest value of e. 

Choosing 9 — tt/2 gives e = 1/4, which is only a mod- 
est improvement on the initial value of e = 1. However 
applying the same network again drives the state closer 
to a pseudo-pure state, and applying the same network r 
times 



tt/2 —) 



I- 
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gives e = l/4 r , so that repeated applications of the net- 
work drive the state asymptotically towards a pseudo- 
pure state. 



r 


£90 


Copt 


Oopt 


1 


0.2500 


0.1437 


77.8° 


2 


0.0625 


0.0294 


99.9° 


3 


0.0156 


0.0008 


96.0° 


4 


0.0039 


4 x 10~ 5 


95.0° 


5 


0.0010 


3 x 10" 6 


94.7° 



TABLE I: The effectiveness of asymptotic preparation of a 
pseudo-pure state in a homonuclear two spin system, com- 
paring optimal choices of the transfer angle 9 with the naive 
choice of 90°. 



For any given initial state it is possible to choose a 
value of 9 which gives more rapid convergence to the 
pseudo-pure state; in general such values must be sought 
by numerical methods. Some results for a homonuclear 
two spin system are given in table U showing that more 
rapid convergence is possible, but these values depend 
both on the choice of initial state and the number of times 
the network is applied, removing the intrinsic robustness 
of the asymptotic approach. 

A further advantage of this approach is that it works 
for any initial state; although the exact results will de- 
pend on the choice of initial state the network will drive 
any state towards the desired pseudo-pure state. Thus 
the same network can be used for homonuclear and het- 
eronuclear spin systems, and for spin systems starting in 
non-thermal states. In heteronuclear systems the state 
will converge more rapidly towards a pseudo-pure state 
if the first gate is applied with the more highly polarised 
spin as the target. For any two spin heteronuclear system 
the optimal angle for the first gate lies within 20° of the 
naive choice of 90°, and so convergence is quite rapid. 

The asymptotic network was implemented using the 
same heteronuclear spin system as before. Fig. [2] shows 
spectra with between r = and r = 4 rounds of the 
preparation sequence Eq. I10[ showing that the state con- 
verges to a pseudo-pure state around r = 3. These spec- 
tra were acquired with a total delay of 18 s between indi- 
vidual scans, only slightly longer than the Ti time of the 
slowly relaxing 13 C spin, so the initial state for these ex- 
periments was not the thermal equilibrium state. As ex- 
pected the asymptotic network is robust against changes 
in the initial state. 



V. EFFECTS OF DECOHERENCE AND 
RELAXATION 

So far we have assumed that the controlled-transfer 
gates will be implemented perfectly, which is unlikely to 
be the case in a real system. Systematic errors can be 
largely suppressed using composite pulses [23 or opti- 
mal control techniques such as GRAPE [llj [2(|, but de- 
coherence and other forms of relaxation will provide a 
fundamental limit. In the context of NMR it is useful to 
distinguish between T2 processes, which cause coherent 
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FIG. 2: Experimental demonstration of asymptotic preparation of pseudo-pure states. Spectra are shown for between r = 
and r = 4 rounds of the preparation sequence Eq. 1101 In this case the spectra were acquired with an delay between experiments 
too short to allow complete relaxation, so the case r = is not simply the thermal spectrum, though it is broadly similar in 
form. Pseudo-pure states prepared with this approach are essentially perfect for r > 3. 



superpositions to decohere, and Ti relaxation which acts 
to restore spin state populations to their equilibrium val- 
ues. Typically T2 is significantly shorter than Ti, and so 
might be expected to dominate, but T2 processes during 
a controlled-transfer gate will only affect the target spin, 
while Ti processes will affect all the spins in the system. 

The effect of T x processes in returning the spin state to 
the equilibrium populations can be modelled using gener- 
alized amplitude damping [24j . and it would in principle 
be possible to largely correct for these effects by sub- 
tly changing the angles in the controlled-transfer gates 
to compensate. A simpler approach, however, is to use 
asymptotic sequences: as these drive the spin system to- 
wards the desired state, whatever the initial state, they 
will automatically compensate for small amounts of re- 
laxation. This analysis ignores the contribution of re- 
laxation to the decoherence of off-diagonal elements dur- 
ing the implementation of a controlled-transfer gate, but 
these effects are indistinguishable from the pure decoher- 
ence effects considered below. 

T2 processes can be modelled in a similar way, and 
their effect is simply to drive the spin system towards 
the maximally mixed state. As such they are mostly 
visible as a reduction in signal size, rather than errors in 
the state prepared. This may be a problem, particularly 
in longer sequences such as those used with asymptotic 
approaches. In our experiments, however, the effects of 
both Ti and T2 processes are small, as the relatively 
large J-coupling leads to comparatively short gate times. 



VI. LARGER SPIN SYSTEMS 

Equivalent networks can be found in larger systems; 
in this case we will only consider homonuclear spin sys- 
tems. With three qubits there are six different controlled- 
transfer gates, leading to a much larger number of pos- 
sible networks. For a network of p gates there are 6 P 
possible networks, each member of which has p different 
transfer angles, 8j. This large group can be immediately 
cut down in three ways. Firstly the networks can be di- 
vided up into groups of six, with the members of each 
group differing only in the labelling of the qubits. This 
is easily handled in a three qubit system by requiring the 
first gate to be controlled by qubit one and to target qubit 
two, implicitly labelling all three qubits. Secondly, there 
is no point in applying the same gate twice in succession, 
as such gate pairs can always be replaced by a single gate 
with a different angle, and so there are only five sensible 
choices for each gate after the first one. Finally, to end 
in a pseudo-pure state, with all populations other than 
the ground state equal, the last gate must normally have 
a transfer angle of 8 P = ir/2. This requirement can only 
be avoided if the network produces a pseudo-pure state 
after the first p—1 gates; in this case, since pseudo-pure 
states are invariant under all controlled-transfer gates, 
the final gate can have any transfer angle. Even in this 
case, however, the network will work with 9 p = tt/2, so 
it is reasonable to impose this as a requirement. 

A numerical search can then be used to minimise the 
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error term e for each of the 5 P_1 networks (five choices for 
each gate except the first) over the p — 1 variable transfer 
angles (the angles of all gates except the last). A brute 
force search is perfectly practical for small values of p. As 
before we assume that the spin system is homonuclear. 
The first perfect solutions (e = within rounding errors) 
are found for networks of five gates, where 104 of the 625 
distinct networks permit essentially perfect pseudo-pure 
states to be prepared. With networks of only four gates, 
the best that can be achieved is e 0.0391, while with 
networks of six gates 1268 of the 3125 distinct networks, 
including all 120 networks containing one copy of each 
of the six gates, permit perfect pseudo-pure states to be 
prepared. 

With systems of three or more qubits it is likely that it 
will not be equally easy to implement controlled-transfer 
gates between all pairs of qubits. An important limiting 
case is when the qubits form a linear chain, with each 
qubit only having couplings to its immediate neighbours, 
so that nearest-neighbour gates are simpler to implement 
than long range gates. This is not a fundamental prob- 
lem, as indirectly coupled gates [25] can be used to im- 
plement long range gates with only moderate overhead, 
but it is useful to investigate whether pseudo-pure states 
can be prepared using networks containing only nearest- 
neighbour controlled-transfer gates. 

For a three qubit system, with four possible nearest- 
neighbour gates, it is once again relatively simple to in- 
vestigate this numerically. The imposition of a chain 
structure means that it is no longer possible to treat 
the three qubits as identical, but a saving of a factor 
of two can be made by noting the symmetry between the 
two possible orientations of the chain (qubit 1 first, and 
qubit 3 first), so in this case there are 2 x 3 P_1 basic net- 
works. Once again the first solutions are found for the 
case p = 5, where two of the 162 distinct networks permit 
perfect pseudo-pure states to be prepared. The explicit 
networks are 
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FIG. 3: The performance of asymptotic networks in a three 
qubit system. The error term e is shown for each of the 120 
distinct networks of six different gates applied to a homonu- 
clear system. The networks are plotted ordered from lowest 
to highest error, and the five "lines" correspond to applying 
the networks from one to five times. 



particular, as previously noted, it is not possible to pre- 
pare a pseudo-pure state of a three qubit system us- 
ing only four controlled-transfer gates even if long-range 
gates are permitted. 



VII. ASYMPTOTIC APPROACHES 

The asymptotic approach can also be used in these 
larger spin systems. We begin by considering the 120 
distinct networks containing one copy of each of the 
six gates, where each gate is assumed to use the same 
transfer-angle, 9 = tt/2. The majority of these networks, 
68, give the same error, e rs 0.1840, but 40 networks give 
lower errors than this, and three achieve the lowest value 
of e « 0.0605, while 12 networks give higher errors, with 
the worst case error, e s» 0.2332, occurring for a single 
network. As usual these results can be improved by ap- 
plying the network repeatedly, and for the main group of 
68 networks the error is given by 



and 



e(r) 



(13) 



tt/2 



(12) 



tt/2 — 



with the angles 9\ = arccos(— 1/7) « 98.2° and 9 2 = 
arccos(— 5/7) « 135.6°. These solutions (together with 
the equivalent pair obtained by swapping qubits one and 
three) appear to be the simplest networks possible; in 



The behaviour of the other networks is more complicated, 
and is shown in outline in Fig. [3] As before, the behaviour 
of these networks can be improved by allowing 8 to vary. 
However initial studies suggest that, unlike the two-qubit 
case, only modest improvements are possible, and so we 
do not pursue this point further. 

Next we turn to the use of asymptotic networks in lin- 
ear chains where only nearest-neighbour gates are per- 
mitted, beginning with the twelve distinct networks of 
four gates. None of these are effective at forming a 
pseudo-pure state, either for the case 9 — tt/2 or when 9 
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is optimised. Applying the networks repeatedly helps a 
little, but is not as effective as was seen in the two-qubit 
case. A likely explanation for this can be seen by exam- 
ining the networks (fTTj) and (|12p : both of these contain 
controllcd-transfer gates with angles of n, which act to 
invert population differences within the system, rather 
than average them out, and this inversion cannot be well 
approximated by any set of averaging gates. 

This can be addressed by relaxing the conditions 
slightly, permitting the use of two angles, namely ir/2 
and 7T, and searching over the 12 distinct patterns of four 
gates with 16 different sets of transfer angles, giving 192 
possible networks. The optimal sequence was found to 
be 



tt/2 



tt/2 



tt/2 



(14) 



for which the error term is given by equation 1131 thus 
this four gate sequence can do as well as the majority of 
six gate asymptotic networks using only angles of 7r/2, 
but it does not perform as well as the best ones. 

Finally we consider the case of networks made up solely 
of nearest neighbour gates with transfer angles of tt/2 and 
7r, but relax the requirement for repeated application of 
a simple pattern of gates, seeking instead the optimal 
network with some total number of gates, p. As before 



such networks occur in mirror symmetric pairs, should 
not contain repeated gates, and should end with a gate 
with a transfer angle of tt/2, resulting in a final list of 
4 x 6 P_2 x 3 distinct useful networks, and we have explored 
these numerically for values of p between four and eight. 
As expected increasing p allows better networks to be 
located, but the best networks for p — 4 and p — 8 are 
no better than the asymptotic networks found previously. 



VIII. CONCLUSIONS 

Controlled-transfer gates provide a simple and effec- 
tive way of preparing pseudo-pure states in small spin 
systems, and experimental implementations confirm that 
as expected these networks give a larger signal intensity 
than conventional approaches. While optimised networks 
provide the simplest approach when the initial state is 
well understood the less-efficient asymptotic approach is 
more robust. In three-spin systems it is possible to find 
both optimised and reasonable asymptotic networks if all 
gates are allowed, but if the spin system is assumed to 
have a chain topology only the optimised approach leads 
to reasonable networks. Initial studies (data not shown) 
suggest that similar results will be obtained in larger spin 
systems. 

Acknowledgments 

We thank the UK EPSRC and BBSRC for financial 
support. 



[1] J. A. Jones, Prog. NMR Spectrosc. 38, 325 (2001). 

[2] D. G. Cory, A. F. Fahmy, and T. F. Havel, Proc. Natl. 

Acad. Sci. USA 94, 1634 (1997). 
[3] D. G. Cory, M. D. Price, and T. F. Havel, Physica D 

120, 82 (1998). 
[4] N. A. Gershenfeld and I. L. Chuang, Science 275, 350 

(1997). 

[5] E. Knill, I. Chuang, and R. Laflamme, Phys. Rev. A 57, 

3348 (1998). 
[6] W. S. Warren, Science 277, 1688 (1997). 
[7] R. R. Ernst, G. Bodenhausen, and A. Wokaun (1987). 
[8] M. Kawamura, R. Sawae, T. Kumaya, K. Takarabe, 

Y. Manmoto, and T. Sakata, Int. J. Quant. Chem. 100, 

1033 (2004). 

[9] M. S. Anwar, J. A. Jones, and S. B. Duckett, Phys. Rev. 

A 73, 022322 (2006). 
[10] X. Peng, X. Zhu, X. Fang, M. Feng, K. Gao, X. Yang, 

and M. Liu, Chem. Phys. Lett. 340, 509 (2001). 
[11] R. Das, T. Mahesh, and A. Kumar, Chem. Phys. Lett. 

369, 8 (2003). 

[12] I. L. Chuang, N. Gershenfeld, M. G. Kubinec, and D. W. 

Leung, Proc. Roy. Soc. Lond. A 454, 447 (1998). 
[13] M. Pravia, E. Fortunato, Y. Weinstein, M. Price, G. Tek- 

lemariam, R. Nelson, Y. Sharf, S. Somaroo, C. Tseng, 

T. Havel, et al., Concept. Magn. Reson. 11, 225 (1999). 
[14] J. Du, T. Durt, P. Zou, H. Li, L. C. Kwek, C. H. Lai, 

C. H. Oh, and A. Ekert, Phys. Rev. Lett. 94, 040505 



(2005). 

[15] J. Zhang, X. Peng, N. Rajendran, and D. Suter, Phys. 

Rev. A 75, 042314 (2007). 
[16] J. A. Jones and E. Knill, J. Magn. Reson. 141, 322 

(1999) . 

[17] E. Knill, R. Laflamme, R. Martinez, and C. Tseng, Na- 
ture 404, 368 (2000). 
[18] M. D. Bowdrey, J. A. Jones, E. Knill, and R. Laflamme, 

Phys. Rev. A 72, 032315 (2005). 
[19] T. Schulte-Herbriiggen, A. Sporl, N. Khaneja, and S. J. 

Glaser, Phys. Rev. A 72, 042331 (2005). 
[20] N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbriiggen, 

and S. J. Glaser, J. Magn. Reson. 172, 296 (2005). 
[21] D. G. Cory, M. D. Price, W. Maas, E. Knill, R. Laflamme, 

W. H. Zurek, T. F. Havel, and S. S. Somaroo, Phys. Rev. 

Lett. 81, 2152 (1998). 
[22] L. Xiao and J. A. Jones, Phys. Rev. A 72, 032326 (2005). 
[23] H. K. Cummins, G. Llewellyn, and J. A. Jones, Phys. 

Rev. A 67, 042308 (2003). 
[24] L. M. K. Vandersypen, M. Steffen, G. Breyta, C. S. Yan- 

noni, M. H. Sherwood, and I. L. Chuang, Nature 414, 

883 (2001). 

[25] D. Collins, K. W. Kim, W. C. Holton, H. Sierzputowska- 
Gracz, and E. O. Stejskal, Phys. Rev. A 62, 022304 

(2000) . 



